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SUMMARY 


Recent developments in supersonic flow theory are applied 
to obtain estimates of the lift-drag ratios that may be 
achieved by aircraft employing swept -back wings. Lift-drag 
ratios greater than 10 to 1 can be maintained up to a Mach 
number of 1.4- by the use of large angles of sweep and high 
aspect ratios. As the speed increases in the supersonic range 
the attainable lift-drag ratios decrease and the g°in due to 
sweepbaclc also appears to diminish. An efficient configuration 
for M= 1.4 would require about o0° sweepback, an aspect 
ratio of 4 and a wing loading of one-third the atmospheric 
pressure. For a wing loading of 50 pounds per square foot the 
cruising altitude would be 60,000 feet and the indicated 
airspeed 290 miles per hour. 


INTRODUCTION 

The work required to propel an airplane a given distance 
in steady flight is equal to its weight times the distance 
travelled divided by the lift-drag ratio of the .airplane. 

Hence the fuel expenditure per mile of flight need not increa.se 
with speed so long as the lift-drag ratio of the airplane can 
be maintained. However, with present shapes a prohibitive 
loss of lift-dr.ag ratio occurs on passing beyond the speed of 
sound and it is evident that a radical change in configuration 
will be necessary for efficient flight at higher speeds. 

The problem of an efficient configuration for flight at 
supersonic speeds was investigated by Busemann in 1935 
(reference 1). Busemann concluded that an improvement in the 
lift-drag ro.tio at supersonic speeds could be obtained by 
sweeping the wing back at an angle just ahead of the Mach 
cone, but failed to recognize the relatively much greo.ter 
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efficiencies obtainable when the wing is swept back behind 
the Kach cone. The change in the type of flow when the wing 
lies inside the Mach cone, and the resulting increase in 
efficiency have been brought out in reference 2. However, 
both reference 1 and. reference 2 are restricted to considera- 
tions of two-dimensional flow and hence aspect— ratio effects 
could not be determined. Recent developments in aerodynamic 
theory have overcome this difficulty, making it possible to 
estimate the lift-drag ra.tio obtainable with practical 
conf igurat ions . 

The present report applies these new theoretical results 
to obtain estimates of the lift-drag ratios that may be 
achieved with an efficient aircraft at supersonic speeds. 

The estimates are all ba.sed on the theory of small disturb- 
ances, first because this is thd only adequate theory 
available , and second because it is reasoned that an aircraft 
producing a large disturbance in the external flow would be 
inherently inefficient. 

At very high Mach numbers even thin bodies and small 
angles of attack cause relatively large pressure disturb- 
ances and consequent heating of the fluid. Here the heating 
effect of friction becomes no longer negligible. (Such 
conditions are likely to be encountered by rockets; however, 
in these cases the efficiency of steady flight may not be of 
primary concern.) The present analysis is therefore limited 
to more moderate speeds where the efficiency in steady flight 
is of primary importance and where it is evident that" such 
efficiency can be achieved by known means. 


FUNDAMENTAL RELATIONS FOR MIN'S LOADING, 

ALTITUDE AND MAXIMUM LIFT-DRAG PATIO 

The lift-drag rs.tio of a conventional airplane depends 
primarily on its external configuration and on the angle of 
attack and 'does not vary greatly with speed provided the 
correct relation between wing loading and altitude is main- 
tained. For maximum efficiency the airplane should be 
flown at that lift coefficient CLont * or which lift-drag 
ratio is a maximum. An increase in" speed, then, necessitates 
an incres.se in altitude, since with fixed lift coefficient 
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(For a complete list of symbols see appendix.) 



( 1 ) 


where the subscript o refers to conditions o/t sea level. 

With lift-drag ratio fixed, higher speed does not involve 
any increase in the thrust required for level flight; this 
thrust is simply 


m 



max 


(2) 


If the propulsive efficiency of the engine docs not drop off 
with altitude, the increase in speed will thus be accomplished 
without any increase in the fuel consumption per mile of 
flight. Furthermore, the increase in speed is not accompanied 
by any significant change in the air loads or pressures on the 
airplane and hence no increase in structural stiffness is 
required. An obvious advantage of this method of increasing 
the cruising: sioeed is that it does not interfere with the 
ability of the airplane to slow down at lower altitudes and 
land on short runways. A more complete discussion of these 
factors will be found in reference 3. 

The altitude and speed of the airplane, of course, cannot 
be increased indefinitely at constant thrust, since eventually 
a critical Mach number will bo exceeded and the lift-drag 
ratio of the airplane will begin to decrease. The limiting 
speed and the corresponding altitude may be determined from 
the relations 

V = (3) 

and 

= CL opt (M 

where Mi is the Mach number at which, for CL = CLopt > 

the drag begins to rise .abruptly, a is the velocity of sound 

and V/S is the wing loading. 
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Equations (3) and (4-) may be combined in the form 



(5) 


where 

p atmospheric pressure at altitude 
7 ratio of specific heats for air (1.4). 

Equation ( 5 ) gives the relation between wing loading and 
atmospheric pressure for maximum speed without loss of aero- 
dynamic efficiency. This condition can hardly be attained, a.t 
low altitudes since with an atmospheric pressure of 2000 pounds 
per square foot, for lu = 0.75 and. the usual values of CL 0P t 
the wing loading required would be of the order of 4-00 pounds^ 
per square foot. At oO, 000 feet, however, the required wing 
loading works out to be the more practical value of 30 pounds 
per square foot. 

Later calculations will show that similar considerations 
apply to supersonic aircraft; that is, the best lift-drag 
ratios are obtained when the wing loading is an appreciable 
fraction of the atmospheric pressure. 

At subsonic speeds it is customary to divide the drag into 
two parts, one the result of friction (including the friction 
dreg of the fuselage) and the other - the induced, drag - the 
result of the lift. The friction drag is considered nearly 
independent of the angle of attack.’ Thus 


where Cd 0 is the drag at zero lift, and for subsonic flow 
equals Cpf, the friction drag. If the velocity and pressure 
disturbances produced, by the airplane are sma.ll, the drag 
arising from the lift will be satisfactorily represented by 
the well-known formula 


CD = Cd 0 + CDi 


( 6 ) 



( 7 ) 
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or 


C I>1 _1_ 

Cl s ttA 


( 3 ) 


and the friction drag will be nearly independent of the angle 
of attach. The lift-drag ratio at any angle is then 


L _ C L 
D CDo+CBi 


( 9 ) 


Solving for the lift coefficient at maximum lift— dra°' ratio 
results in 0 



( 10 ) 


and therefore 


f-' = I /~ 1 

l^/'max 2 / CDo (CDg/CL 2 ) 

V 


( 11 ) 


In calculating lift-drag ratio for supersonic speed the 
drag may again be divided into' two components, one inde- 
pendent of the lift and one proportional to the square of the 
Iit"t coefficient, Hoiirevcr, in this case the drop* at zero lift 
includes a pressure drag which varies with the thickness of 
the body or wing. Also, at supersonic speeds, the drav duo 
to lift can no longer properly be called" "induced drae." At 
subsonic speeds the drag arising from the lift can be traced 
to the influence of the trailing vertex wake on the wine: - 
hence the designation "induced." At supersonic speeds, ~ 
however, the forward influence of the wake usually constitutes 
only a small part of the drag arising from the lift and hence 
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the term "induced drag" does not seem appropriate. Different 
divisions of the drag due to lift into components of wave 
drag and induced drag have been proposed, but the proportions 
allotted in any particular case depend on the method of calcu- 
lation employed. In the present report the drag is calculated 
by integrating the pressure distribution in the neighborhood 
of the body and in this ca3e Cpj_ appears simply as a 

pressure drag proportional to the square of the lift coeffi- 
cient. The subscript i is retained to identify the law of 
variation with that of the induced drag at subsonic speeds. 

Then, for comparison with the subsonic case, we may write 


CD 0 = Cp f + Cp t 


and 


Cp 


•i 


/'CDin 

VcITV 


cl 


2 


where Cpf is the total friction drag, and CDt the total 
thickness drag, due to wing and fuselage. The factor 
CDj_/CL 2 bears no simple relation to the aspect ratio as it 
docs in the subsonic case, but is a complex function of the 
wing plan form and load distribution. 

With the values of C.p 0 and Cpp/Ci/, revised for super- 
sonic conditions, eouations (10) and Til) for the optimum lift 
coefficient and maximum value of the lift-drag ratio remain 
valid. Maximum L/D is obtained when the drag due to lift is 
equal to the drag at zero lift. 

DRAG AT ZERO LIFT 
Thickness Drag of Wings 

The thickness drag of the wing may bo calculated, by the 
methods of reference 4 or 5* Fig'ure 1 shows the variation of 
thickness drag with Mach number calculated by the method of 
reference 5 for a rectangular wing and for several swept- 
back wings. In these cases the e.irfoil is of symmetrical 
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biconvex section 5 percent thick, 
swept-back airfoils were obtained 
curve for the rectangular airfoil 


“7 

( 


I 


The results for the 
from reference 6. The 

n r> »-• 


, _ is the same- as that given 

for the infinite wing by the Ackcret theory since, as has 
been demonstrated oy J. N . Nielsen in an unpublished ampli- 
cation of the same method, the integrated effect -of removing 
outboard portions of the wing on the "drag’ of the remainder of 
tn.u wing is zero, at least so Iona as the Mach cone from each 
oip does not intersect the opposite tin. The- thickness dram 
coefficient of the rectangular airfoil is therefore- "constant 
when 


A Vms_T> i 


where 


A aspect ratio 

When the wing is swept well behind the Kach cone the 
ilow over most of the wing is of the subsonic, tvpe. (See 
reference 2. ) Jhe pressure drag is small and may be 
-treated to departures from the subsonic type of flow in 
the region of the root section. In this condition the' 
outboara sections of the wing have little or no dra- r. nc \ 

esS?ct t ratio aS At°ff i ’ CiCnt is , inverse ly proportional' to the 
aspect ratio. At higner speeds, when the Mach male 

approaches the leading-edge angle, the distribution of drag 
changes and tne drag coefficient increases rapidly -orrtic- 
ularly on the outboard sections. If the leading edge is 
too near tne Kach cone, the drag of the swept wine will 
exceed, that of the straight wing. 

. * igure 2 shows a plot similar' to figure 1 of the 
variation of drag with Mach number for tapered swept-back 
airfoils. Those results were obtained by X, HarmoHs of 
Langley memorial Aeronautical Laboratory usinv the method 
o: reference 5. An extensive series of' calculations for ' 

Win ff na f been ^en recently by Stewart and Puckett 
(. e.erence 7). In order to simplify the calculations a 


double-wedge section was assumed, though it 


lB not to be 


Data 


on file at Ames Laborato - *’*'. 
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supposed that such a section would be desirable in. practice. 
Here the angle of sweepback (60°) is that of the midchOi 1 
line of the~airf oils , which is a.lso the line of maximum thick- 
ness. The sharp rises in drp.g coefficient near M = l.p- and 
M = I .71 occur" when the Mach angle approaches the angle of 
the tre.ilinec edge. Evidently s.ll generators of ohe win*, 
surfa.ee must lie behind the Mach lines to insure fo.vorablu 
drp.g values. 


According to the thin airfoil theory the calculated flow 
for a. given airfoil plan form and Mach number will actually 
be similar to the flow over another plan form at a different 
Mach number, provided the two plan forms are oriented similarly 
with respect to the corresponding Mach lines. This relation 
may be preserved by changing the x coordinates of the pla.n 
form (fig. 3 ) in the proportion that the x coordinates of 
the Mach lines are changed, that is, as —1 • i 1 or plan forms 
having similar flow patterns the ratio 

m _ cot leading-edge angle . . wlll be constant, 
cot sweep angle of Mach lines 

The aspect ratio will then vary with Mach number according to 


A -1 = constant (12) 

s.nd thin r.irfoil theory shows the.t the drag coefficient will 
be proportiona .1 to 

(t/c ) 2 

V i: 2 -i" 


or 


ODt VT^l ~ (t/c ) 2 


(13) 


where t/c is the thickness- chord re.tio mee.sured in the stream 
direction. Figure 4 - shows the coefficient 


CDt */&=! 
(t/cT 2 
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plotted against 


m 


for the constant-chord "biconvex airfoils. 


Friction Drag of ’Ting 

In general, the friction drag of the wing will he of the 
sane order of magnitude as the thickness drag. At. very high 
speeds a. considerable amount of heat is generated in the 
boundary layer and the resultant temperature variation affects 
the magnitude of the skin friction. For moderate supersonic 
speeds,' however the heating effect is not large a.nd the normal 
relation of shin friction to Reynolds number will not be 
greatly modified. 

For present purposes a conservative value of = *“** ^06 

corresponding to a turbulent boundary layer at a Reynolds 
number of 10 7 has been used. 


Drag of Fuselage 

A method for calculating the wave drag of a ^slender 
fuselage at supersonic speeds was given by von Karman in 1935 
(reference gj ." This method was applied in reference 9 to a 
series of bodies of parabolic arc shape and estimates of the 
friction drag added to obtain total drag. More recently the 
calculations' of Haach (reference 10), Sears (reference 11), 
and Lighthill (reference ,12) have become available. These 
investigators apply Kerman* s method to the determination of 
body forms having a minimum wave drag for certain conditions. 
The minimum problem is solved for three cases: viz, T, given 
volume and given lsngtlij II, given Length a.nd given dia.iue oer , 
and III, given diameter and given volume. 

The following equations may "be obtained 2 from Haach’s. 
report (reference 10). The length l is so chosen tnat tne 
body lies between +1 and -1 on the x-axis; r/r 0 is the 
radius at anv station in terms of the maximum radius r 0 and 


a is the fronta.l area rrr o' 


;he volume is given in terms 


of the volume of the circumscribed cylinder, and the crag 
coefficient, which does not include friction, is given in 


*The formulas given as the final relations in the 
in error . However, the correct relations ca.n eau 

derived from the preceding equations. 


report 
ily be. 


are 
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terms of the frontal area, 
ratio, diameter /length. 


The fa.ctjor d/ 1 is the 


Case I: 


Given length, 


given voluir 


e 



Volume 



rr 0 


2 


c 




Case II: Given length, given diameter 


( — — ) 2 = VI- x 2 - x 2 cosh -1 — 
\ r 0 J x 

Volume s I l nr 2 

6 



Case III: Given diameter, given volume 

(— S ) Z = 3V1-X 2 - 2 ) 3 - 3x 2 cosh- 1 — 

\roJ x 

Volume = 5- l nr 2 
5 ° 

3 _ 2 ( £\ 2 
° D = 2 11 


fineness 


s (m-) 


> ( 15 ) 


) (16) 

! 

i 

i 

) 


Figure 5 shows the body sha.pes computed from these formulas. 
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Although the waTs drag diminishes with increasing 1 slen- 
derness, the friction clrae for a given volume or a given 
cross section tends tc increase because of the greater 
surface area. With usual values of the friction coefficient 
o. favorable balance between the two components rccuires such 
a. slender body that in most cases the dimensions will actually 
be .governed by the minimum allowable cross section. For a 
slender body the surface area and hence the friction drag 
associated with a, given cross section is proportional to l/d. 


It follows 


while the wave drag is proportions.! to (d/l) 
that the total drag will be a minimum when the slenderness 
ra.tio is such that the friction dram is twice the 


ave drag. 


It will be noted that the body shape 
has very little more drag than the case I 


for case I actually 
body of the same 


diameter, and since body I has a. grer 
seems a. logical choice for pre.ctical 
the wave drag for ca.se I and also th- 
skin- friction coefficient of C.0021, 


ter useful volume, it 
design. Figure C shows 
total drag, based on a. 
as a function of the 


fineness ratio. The value of this friction coef 
obtained from reference 13, and corresponds to a 
lent boundary layer and a Reynolds number of 10 3 
friction coefficient the optimum fineness ratio 


ficiont we. 8 
fully turbu- 
. With this 
is s/bout 16 


to 1. 


DRAG 


DUE 


LIFT 


The dra.g due to lift is estimated, from theoretical 
solutions for the supersonic flow over thin lifting surfaces. 
Theoretical solutions ere known for cases in which the 
lifting surface is curved a.nd twisted in such a. way a.s to 
support a. uniform load (reference lL) and, for certain 
rectangular, triangular, or taoored flat surfaces (references 
14 - and "l5 ) . 


Uniformly Loaded. Surfa.co 


The solution for the 
derived by methods similar 
for the nonlifting airfoil 
to thickness on an airfoil 
by superposing the effects 
acceleration potential fie 


uniformly loaded, surface ma.y be 
to those described in reference 5 
. In that report the pressure due 
obi i era e to the stream was obtained 
of oblique line sources in the 
Id. The effect of a line source Is 
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to cause a deflection of the stream linos- crossing the source 
like the deflection caused by a thin wedge-shaped body; ’that 
is, the line source is followed by an area over which • the 
vertical velocity w is constant and of opposite sign above 
and below the chord plane. 

Similarly, an obliouc vortex gives rise to. a. constant 
difference in the horizontal velocity increment u, and there 
fore in the pressure, above and below the plane of flow cross- 
ing the vortex. The corresponding w for a semi-infinite 
vortex is given by 


w = — .■Vl-m cosh -1 — ~ - cosh -1 -3L~\ ( 17 ) 

rrra ^ . y'i jy \) 


where x 1 = x - my and y' denotes the absolute value of 
y - mx (The geometry of the figure has been adjusted, as 
described in the preceding section and reference 5, to 
correspond to the cs.se in which the Hack a.nscle is -S°\ that is 

= i. ) 


The shape of the surface and the 
related to the velocity increments by 


constant pressure are 
the following formulas: 


c.z _ w 
dx V 


and 


( 13 ) 


Ap = 2u 
q V 


( 19 ) 


Thus the camber of a triangular airfoil shaped to support a 
uniform load (fig. 7) may be obtained by superimposing two 
oblique vortic-os to form a V coinciding with the leading 
edge of the triangle. Integration of equation ( 17 ) for this 
case yields: 
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z = 


Cl f v / I^¥ 5 ' / - 


* 3C * tr I — ~ 1_ 


s-TTm 


/ y 1 cosil — -y* cosh 

m v, |y«j 


— ) 

|y’l J 


\ (20) 


-2 (x cosh 


- 1 x 


1 

\ 


-'•/ x 2 - y ^ ) j 


where 


x 1 = x + my and y 1 = y + mx 


To obtain a, lifting 1 surface of finite chord it is necessary 
to introduce a negative V-shaped vortex at the desired, chord, 
length downstream. (See fig. S. ) Through the use of a finite 
number of straight vertex segments any ulan form bounded by 
straight lines can be obtained. 

The variation of w over the area enclosed by the 
vortex segments not only gives the camber and twist of the 
surface required, to support a uniform load, but also can be 
used to calculate the drag arising from the lift. It can be 
seen that, since the pressure distribution is uniform over 
the section, the resultant force will lie in a direction at 
right angles to the chord line, or the line .joining the 
leading and. trailing edge, regardless of the camber of the 
surface. Hence the angle of attack of the chord line at any 
section times -the lift gives the drag d.ue to lift at that 
section. 

In case the leading edge of the airfoil is ahead of the 
Mach cone the uniformly loaded surface is flat over portions 
of the wing not influenced b3 r the root or the tip, as is given 
by the Ackeret theory. More interesting cases are those in 
which the leading edge's are swept behind the Mach cone. 

In the case of the swept-ba.ck wing it is found, that the 
angle of attack has a logarithmic infinity at the center 
section. Hence the wing would require an infinite twist to 
maintain the uniform load across this section. At a distance 
froin the center section the shape of the lifting surface 
resembles that of the familiar "constant load mean line" used 
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for subsonic airfoils. The twist and hence the section drag 
disappear rapidly with distance from the center section. 

There is consequently a narked reduction of drag coefficient 
with increasing aspect ratio, just as in the case of the drag 
due to thickness. 

The infinite twist required at the root section, of 
course, makes the construction of such a wing impractical. 

We may conclude that in a practical wing there will be some 
falling off of the lift ©.cross the center section, and calcu- 
lations of the lift distribution for flat surfaces show such 
a loss. The uniformly loaded airfoil gives a useful picture 
of the variation of drag with plan form, however. In spite of 
the fact that the local drag coefficient at the root section 
tends toward infinity, the integrated or over-all drag 
coefficient of the swept-back wing is finite and at reasonable 
aspect ratios is considerably lower than that of the flat 
unswept wing. 

Figure 9 shows the coefficient of drag due to lift 
Cpi/Ci, ""for a series of uniformly loaded airfoils having a 
constant chord and varying degrees of sweep. To simplify the 
calculations, an approximation was made for the effect of the 
wing tip. With the’ tip cut off parallel to the direction of 
flight a large twist Tr ould theoretically have been required to 
maintain the uniform load right out to the tip. Instead of 
calculating this additional twist at the tip, the shape of the 
infinite wing with uniform load was assumed without modifica- 
tion and a loss in lift within the Ilach cone originating at 
ea.ch tip was taken into account. Since the lift will have the 
full value along the boundary of the cone and will fall to 
zero at the tip, an average value of half the full load was. 
used over this region. Since the effect of this approximation 
to a tip effect on the total drag value was small, any error 
involved in the approximation must also be small. If the tip 
were cut off along the Mach lines, slightly lo^er values of 
Cpi/Cb would have been obtained. 

Figure 9 shows that the values of Cp^/Cp, super- 
sonic speed are in general higher than the^value corresponding 
to the same aspect ratio at subsonic speed but approach this 
value as the angle of sweep is increased 'i.e., as m->0). 

The practical difficulty of maintaining a given aspect ratio 
of course increases as the angle of sweep is increased. 
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Flat Lifting Surfaces 

Rectangular plan form.- For a flat rectangular wing of 
infinite aspect ratio the Ackeret theory gives 


6C L _4 

6 a 


(21) 


Since the lift is at right angles to the chord, 


= Ct. x a 


( 22 ) 


and 


1 °Di 

C L - 


0.25 


(23) 


At a Mach number of 1.4- this value is nearly five tines the 
drag due to lift of a subsonic airfoil of aspect ratio 6. 

If the rr ing has a finite aspect ratio there will be a 
reduction of lift at. the tip and a consequent reduction in 
6CL/6a from the value given by equation (21). The distri- 
bution of lift over the tip of a flat rectangular wing has 
been calculated by Busemann (reference l4) . The lift over the 
portion of the rr ing between the tip Mach cones (fig. 10) is 
constant and equal to that given by the Ackeret theory. 

Within either tip cone the lift pressure falls from this value 
to zero at the tip. If y/x represents the fractional 
distance from the* tip toward the Mach line at a given chord- 
wise position, then the lift pressure varies according to 
the function 


Jh_ 


cos 


— 1 


0 



(24) 


l6 
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Superposing the effects of the two tip cones where they. • 
overlap and integrating the pressure over the whole wing gives 
for the lift coefficient 


Cl 


4a AVh 2 -l -4 
Vh 2 -! A Vh 2 - 1 


(25) 


and for the drag due to lift 


VTPCT 


'C L 2 - 


) = 


Aa/T^NT 

-2 


(26) 


Busemann 1 s solution is valid for A vM 2 -l > 1.0, that is, 
so long as the Mach cone from one tip does ""not cross o ver the 
opposite tip. It is interesting to note that ^hen AVM 2 ~1=1.0 
the lift falls to zero along the Tr hole trailing edge and the 
span load distribution is elliptical, as shown in reference l6 
for airfoils of very low aspect ratio. 

Triangular plan form.- Formulas for the lift distribution 
and bCL/da for a flat triangular airfoil behind the Mach 
cone have been given recently by Stewart (reference 15)* 
Stewart finds that the lift distribution as predicted from 
elementary considerations for very slender triangles (refer- 
ence 16 ) actually holds for all leading-edge angles until the 
leading ed.ge touches the Mach cone. Stewart finds also 


VPW ib = 2m (27) 

oa E 


where E = E (Jl-n y ) is the elliptic integral. 

In the case of the flat surface with the leading edge 
behind the Mach cone, the chordwise lift distribution has an 
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infinite value at the leading edge just as it does in the 
subsonic case. Here the resultant force will be inclined 
forward rela.tive to the chord plane because of the suction 
force at the leading edge. 

The drag due to lift for the flat triangular airfoil 
was evalus.ted by setting up the complex expression for the 
velocity field u and. w by means of Busemann 1 s method 
(reference 1^). The drag was then calcule.ted from the 
formula 



(2d) 


by integrating around a contour c a short distance away 
from the airfoil surface and enclosing the singularity at 
the leading edge. The result is 


1 2E - */l - m 2 

Cl7 K - mn 


( 29 ) 


A similar formula, has been given recently bv P. Hayes 
(reference 17 ) . 

In this formula, the first term represents a drag equal 
to the lift times the angle of attack, a.nd the second, term 
represents the thrust at the leading edge. It is noted that 
this latter term disappears progressively as the edge 
approaches the Mach cone ( i . e . , 'Avto s -1— >4 ) . At the 'other 
limit, the slender triangle near the center of the Mach cone 
E_>1, and ~ ’ 


lg 
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as in reference l6. Although- the theory shows a forward 
thrust on the thin plate with a sharp edge, it is not to be. 
expected that this characteristic will be realized in practice 
unless the leading edge is given a finite radius or caiaber. 

Tapered plan form.— The theoretical lift distnoutxon 
for a flat un tape red swept-back wing with the leading edge 
behind the Mach cone has not yet been determined. However, the 
solution for the flat triangular wing may be readily extended 
to include a special family of tapered wings. This extension 
is based on the fact that an area of the triangular wing may- 
be removed by making cuts alone Mach lines witnout afiocting 
the flow over the area remaining ahead of the cuts. in 
particular, the removal of such area will not affect the 
suction force on the leading edge, as long as the area removed 
does not include any of the leading edge so that the 
coefficient of thrust will be increased as area is cut away . 
Evidently the most efficient members of this family of air- 
foils are those in which the. maximum area is cut out of the 
triangle, that is, the wing is tapered to a point. (See 
fig. II). 

With the trailing edge fixed. at the Mach angle, the angle 
of taper and hence the aspect ..ratio of these wings varies with 
the ancle of sweep in such a manner that 


a = 


lirr 

I-m 


(31) 


as the leading edge approaches the Mach cone m— >1.0 
and the aspect ratio approaches infinity. 


The lift-curve slope of these airfoils is determined 
simply by integrating the pressure distribution for the 
triangular airfoil over the appropriate area, 
tion gives 


for 

The calcula- 


1 / cos 1 (-m) + m ,/ 1-m 
£>o Vl-m 2 E v 1 + m 1+n 




( 32 ) 


md for the drag due to lift 
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_1 ^£i = 1 . _ / 1+m 1 IT 2 

VM 2 - 1 cl 2 Vm 2 -! 6Cl v l-m Wi (33) 

> . 6 a 

where N is the ratio of the lift— curve slope of the 
triangular airfoil to that of the tapered airfoil, that is, 


N ' = . (tt/SV l-n g | ( 3 ![.) 

cos" 1 (-mV 

1+m ' V l + ro 


Equations (31) through (3M-) apply to the case of the 
wing tapered to a point. In subsonic flow such extreme- taper 
is known to lead to high local lift coefficients over' the 
tip portions and to the possibility of tip stalling even- at 
moderate lift coefficients. A similar tendency is evident 
at supersonic speeds; in fact, the section lift coefficients 
tend toward infinity at the pointed tip. Hence the extreme 
taper should not be used in practice and value of Cp^/Cps 
calculated for these cases will be somewhat optimistic. 


Comparison of Lift and Drag. Values for Flat Surfaces 


Curves showing the variation of lift-curve slope with 
Ma.ch number and a.spect ratio for the rectangular, triangular, 
a.nd tapered airfoils are shown in figure 127 • At A/h 2 -! = 4 
the leading edge of the triangular airfoil touches the Mach 
cone and, as shown by Puckett - (reference 4), the lift charac- 
teristics at higher a.spect ratios are identical with those of 
-a. rectangular airfoil of infinite aspect ratio. 


The drag due to lift versus aspect ratio for the various 
flat wings is shown in figure 13. According to the Ackeret 
theory 


1 C Pi _ 1 
JW-T Cl 2 ^ 


(35) 


and it is - to be noted- that both the rectangular a.nd the 
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triangular airfoils approach. this value at higher aspect ratios, 
At 




(36) 


the leading edges of the triangular airfoil are behind, the' h'ach 
cone and the drag due to lift is reduced somewhat because of 
.the suction on the lead.ing edge. However, the really favor- 
able values of Cpj_/CL 2 are obtained, only with the swept-bacii 


wings of relatively high aspect ratio. The fact that the 
values for the flat pointed wings agree with those for the 
cambered, untapered.. 'air foils shown on figure 9 is an indi- 
cation that the drag, due to lift is primarily a function of 
sweepback and aspect ratio. 


RESULTS 

The total drag of the supersonic aircraft can now be 
estimated by adding up the components thus far considered 
with an allowance for the friction drag cf the wing and a. 
small allowance for the tail surfaces. 

Since the lift-drag ratio increases with increasing slen- 
derness of the wing, it is necessary to establish some standard, 
of slenderness to obtain comparative values. A rough measure 
of the structural stiffness of a wing is the maximum spar 
depth at the wing root divided by the distance, measured along 
the spar, to the centroid of area of the wing. A value of 
l/l5 seems to be about the limit of present-day construe tion. 


Airplane with Constant Chord Swept -back Wing 

Figure 14 shows lift-drag, ratios obtainable at VL = 1.4 as 
a function of m with a configuration embodying the constant 
chord, uniform lift airfoil and a type^I body of 15 to 1 fine- 
ness ratio, An allowance of CPf= O.OOo.'was made for the 
friction drag on the wing and a value eoual to 10 percent of the 
wing drag was allotted, to the vertical tail. No horizontal 
tail is shown, since it is not clear that such b. tail would be 
reouired with this configuration. The frontal area of the body 
was assumed, to be 4 percent of the wing area. The drag and lift 
of the wing we re assumed, to carry across the center sections. 
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without being modified by the presence of the body. 

The airfoil shape is obtained by superimposing a Para- 
bolic arc thickness distribution upon -a cambered and tvisted 
surface designed, as discussed earlier, to support a uniform 
load. The varion in sweepback, or m, in this case was 
assumed to be obtained by rotating the wing panels without 
changing their length-width ratio, hence the aspect ratio 
varies with sweep as shown. The wing in the unswept position 
would have an aspect ratio of 12. If the ratio of the root 
thickness to the spar length from the root section to the 
centroid of the wing panel is 1 / 15 , the thickness- chord ratio 
of the unswept wing would 0.2. The same wing, rotated 
through 6o° (m = 0.577) has an aspect ratio of 3 and a t/c 
(c measured parallel to the stream} of 0.1. 

The calculations for the uniformly loaded wings show 
higher lift-drag ratios for still higher aspect ratios and 
greater thickness-chord ratios, but it is doubtful that the 
calculations based on the theory of ema.ll disturbances apply, 
in these cases. 

Because of the higher aspect ratios attainable with a 
tapered wing it is found that these configurations are more 
efficient than the constant chord wings, and therefore they 
will be discussed in somewhat greater detail. 


Airplane with Tapered, "ing 

Figure 15 shows the lift-drag ratios obtainable at H = l.i 
with the flat pointed wing. The proportions of fuselage and 
tail are the same as in the preceding case, and the same value 
of Cp-f was used. The calculations were made assuming double- 
wedge sections (as in fig. 2), but an approximate correction 
factor of i/3 was inserted into the thickness drag to take 
account of the greater average slope of the biconvex profile. 
The quantity i/3 is the ratio of the wave drag of the biconvex 
section to that of the double-wedge section in two-dimensional 
flow. The maximum wing thickness in each case is again 1/15 
the distance from the wing root to the centroid, of area of the 
half wing. The variation of t/c (streamwise) is shown in 
figure 15 . Since the trailing -edge angle is fixed on the Mach 
line, the aspect ratio increases indefinitely as m approaches 
1, according to equation (31). The optimum lift coefficient 
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calculated by equation (10), if? a.lso shown in figure 15 * 

Figure l6 shows the variation of (L/D)m&x x,r it:h n at 
different Mach numbers . It is noted that the optimum value cf m 
is different for different Mach numbers; there appears to be 
no fixed relation between the- sweepback angle and the Mach 
angle. Evidently relatively greater sweepback angles should 
be used at smaller Mach numbers. The optimum values of m will, 
of course, be influenced by the magnitude of the friction drag. 


DISCUSSION CF RESULTS 
Effect of Plan Form 

Figure 17 shows the lift-drag ratios replotted against 
aspect ratio and compared with values estimated for a straight 
wing— body combination. It will be noted, that up to M = 2 the 
swept— back wing is much more efficient than the straight wing. 
The difference is « mailer at the higher Mach numbers, however, 
and the advantages of sweepback at very high Mach numbers may 
be questioned. In each case the efficiency diminishes with Mach 
number. 

Although the configurations shown in figures. lk and 15 
appear from the calculations to give the best lift-drag ratios, 
it is not to be assumed that these configurations are actually 
the most suitable for practical use. In Practice the wing must 
of course have a finite tip- chord and may also require some 
camber or twist to avoid the high concentration of load, near 
the tips. Also, as has been previously remarked, the location 
of the trailing edge on, rather than behind, the Mach lines 
was chiefly a computational device. It is probable that a 
greater sweep of the trailing edge would, be desirable. 

Such modifications will of course cause changes in the 
lift-drag ratio. However, it is believed, that the highest lift- 
drag values shown can actually be approached with practical 
configurations. The theoretical values of Crg /Cl* for the 
wing with its trailing ed.ge along the Mach cone are somewhat 
more . favorable than the values to be expected. T - r ith a. wing 
having its trailing edge behind the Mach cone. On the other 
hand, the location of the trailing edge along the Mach line is 
unfavorable from the standpoint of thickness drag, as shown by 
figure 2. Hence the net effect of tra.il ing-edge location on 
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(L/D)max is not expected to be very pronounced. The benefi- 
cial effects of tapering the wing indicated by figure 17 may 
also be assumed to hold qualitatively for more moderate degree 
of taper. 


Airfoil Section 


No attempt was ms.de in the analysis to find an optimum 
airfoil profile. The section assumed for the calculations 
has a parabolic thickness distribution. In practice, as 
previously mentioned, it would be necessary "to round or 
camber the leading edge to achieve the predicted values of 
Cdi/Cl 2 » It might also be advantageous to use (with the 
tapered plan form) a cusped trailing edge. This device would 
ens.ble the designer to te.ke advantage of the high lift to be 
obtained by pls.cing the trailing edge along the Mach lines, 
while effectively giving the thickness distribution a gres.ter 
angle of sweep and thus averting the Is.rge wave drag which 
arises when the generators of the thickness distribution are 
too near the Mach lines. 


Friction Drag 

The allowance of 0.006 made for the friction drag coeffi- 
cient of the wing corresponds to a turbulent boundary layer 
at a Reynolds number of 10 7 . The assumption of turbulent 
friction for both wing and fuselage is believed to be 
conservative, since there' are indications that large areas of 
laminar flow can be achieved at supersonic speeds." The 
importance of maintaining laminar flow or otherwise reducing 
the friction can be seen from the magnitudes of the various" 
drag components with the best configuration (fig. 15, m = 0.5) 
at a Mach number of 1.4. The various components are shewn 
in the following table: 


(1) 

Thickness drag of wing 

o.oo4i 

(2) 

Friction drag of wing 

.0060 

(3) 

Thickness drag of body 

.0020 

(4-) 

Friction drag of body 

.0036 

( 5 ) 

Drag of vertical tail 

. 0010 

Total 

drag at zero lift : 0p o = 

.0167 


Drag due to lift Cp* = 

.0167 

Total 

drag 

0.0JS4 

Friction drag (2) + (4) 

0.0096 
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Note that the friction drag is more tha,n 5 ^ percent of the 
total drag at zero lift. In this case the maximum lift-drag 
ratio is 10.7,* with completely laminar flow the ratio would 
increo.se to about 15* 

Optimum ■.Ting Loading and Altitude 

The analysis indice.tes that reasonably good aerodynamic 
efficiencies are obtainable up to Mach numbers of 1 . 3 * At 
M = 1.4 the best configuration studied should operate near a 
lift coefficient of 0 . 3 R. From equation ( 5 ), the wing loading 
for this case works out to be about one— ho.lf the atmospheric 
pressure. This pressure disturbance co.n no longer be considered 
small a.nd the question .arises as to whether the linearized 
theory can be considered applicable in this ca.se. No accurate 
analysis of this limitation can be given a.t present. However, 
a.n approximate criterion can be d.ed.uced by comparing the flow 
over the swept-back wing with the two-dimensional subsonic flow 
over a wing section a.t the same component Mach number a.s 
suggested in reference 2 . -.Then this comparison is made for 
configurations near the optimum in figure 15 it is found, that 
the wing sections are operating beyond their critical Mach 
numbers at the ind.ica.ted optimum lift coefficient. Thus it 
appears that the optimum lift coefficient will actually be 
smaller than is indicated by. the linearized theory. For the 
best configuration a.t M = 1 . 4 , it appears that tuc optimum 
lift coefficient may be nearer 0.25 than the 0.357 indicated 
by figure 15. In this case the L/Dmax will be diminished 
from 10.7 to 10, and the optimum wing; loading from one-half to 
approximately one-third atmospheric pressure. At sea. level 
the wing loading required 'would be 700 pounds per square foot, 
but for operation at 60,000 feet the much more reasonable figure 
of RO rounds per sq.ua.re foot is obtained., .At this altitude the 
true airspeed is 900 miles per hour a.nd the indicated airspeed 
290 miles per hour. 


Ames Aeronautical Laboratory , 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., May 19 ^ 7 • 
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APPENDIX 

SYMBOLS 

V flight velocity 

Vo velocity at sea level 

p air density 

Po density at sea level 

p atmospheric pressure 

T thrust 

¥ weight 

L lift 

D drag 

M Mach number 

a velocity of sound 

Cl lift coefficient ( —If ^ 

S wing area 

7 ra.tio of specific heats ( 7 = 1.4 for air) 

Cp draff coefficient ( D '\ 

Vp/2 V 2 S J 

CDo drag coefficient e.t zero lift 

Cpt coefficient of drag due to thickness 

Cp^ coefficient of dra.g due to lift 

Cpf friction drag coefficient 
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aspect ratio 

wing span (perpendicular to direction of flight) 
coordinate along direction of flight 

parameter indicating relative alone of wing leading edge 

/ _ cotangent 3 eadinn-eclge angle \ 

\ m = cotangent' 'sweep ' angle of* 'uach' lines ) 

thickness of wing at midchord 


wing chord 
length of fusela.ge 
radius of fuselage 
maximum radius 

maximum diameter of fuselage (d = 2r 0 ) 
small vertical velocity disturbance 
small horizontal velocity disturbance 
lateral (spam r ise) coordinate 
x - my x' = x + my 


y - mx y 1 1 = y + mx 

vertical coordinate of wing camber line 
angle of attack 

complete elliptic integral of the second kind 

ratio of lift-curve slope of triangula.r airfoil to 
that of the tapered airfoil 
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Fig. 9 




Figure 10.- Lift distribution over flat rectangular surface (ref. 7). 
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Figure 11.- Lift distribution over swept-back tapered wing with trailing edges along 
Mach lines . 


NACA TN No. 1350 Fig. 11 




NACA TN No. 1350 Fig. 12 




F/gure /3 - D/ 'OCJ r/cje fo //Ff~ 09 a funcf/on of 05/^e.cf ratio eo'nyp 


Fig. 13 NACA TN No. 1350 


NACA TN No. 1350 


Fig. 14 
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Fig. 16 
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